Higher twist effects in charmed-strange z/DIS diffraction 
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Abstract 



The non-conservation of charmed-strange current in the neutrino deep inelastic scattering 
r^, ' (z/DIS) strongly affects the longitudinal structure function, Fl, at small values of Bjorken x. 

The corresponding correction to Fl is a higher twist effect enhanced at small- a: by the rapidly 
growing gluon density factor. As a result, the component of Fl induced by the charmed- 
strange current prevails over the light-quark component and dominates Fl = F[ s + F^ d at 
x < 0.01 and Q 2 ~ m 2 c . The color dipole analysis clarifies the physics behind the phenomenon 
and provides a quantitative estimate of the effect. 
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1 Introduction 



Weak currents are not conserved. For the light flavor currents the hypothesis of the par- 
tial conservation of the axial-vector current (PCAC) [1] provides quantitative measure of the 
charged current non-conservation (CCNC) effect. [2]. The non-conservation of the charm and 
strangeness changing (cs) current is not constrained by PCAC. Here we focus on manifesta- 
tions of the cs current non- conservation in small-x neutrino DIS. At small x the color dipole 
(CD) approach to QCD [3, 4] proved to be very effective. Within this approach it is natural 
to quantify the effect of CCNC in terms of the light cone wave functions (LCWF) 1 

M> ~ ge v f/AE, (1) 

where j v = c(k)^ u (l — 7s)s(p), AE = E q — E p — E k and e v is the four-vector of the so-called 
longitudinal polarization of the W^-boson with the four-momentum q. Notice that e v — > q u /Q 
for Q 2 = -q 2 -> 0. 

The observable which is highly sensitive to the CCNC effects is the longitudinal structure 
function Fl(x, Q 2 ) related, within the CD approach , to the quantum mechanical expectation 
value of the color dipole cross section, 

F L ~ Q 2 (*|<t|*). (2) 

Our finding is that the higher twist correction to F L arising from the cs current non- conservation 
appears to be enhanced at small x by the BFKL [6] gluon density factor, 

*r~4(±) A . (3) 



Q 2 \x 

The color dipole analysis reveals mechanism of enhancement: the ordering of dipole sizes 



(m 2 c + Q 2 )' L <r 2 <m- 2 
typical of the Double Leading Log Approximation (DLLA) and the multiplication of log's like 

a s \og[(m 2 c + Q 2 )/^ G )\og(l/x) 



1 Preliminary results have been reported at the Diffraction'08 Workshop [5] 
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to higher orders of perturbative QCD. As a result, the component Fff induced by the charmed- 
strange current 

F L = F u L d + F C L S (4) 
grows rapidly to small- x and dominates F L at Q 2 < m 2 [7, 8] . 

2 CCNC in terms of LCWF 

In the CD approach to small- a: z/DIS [9] the responsibility for the quark current non-conservation 
takes the light-cone wave function of the quark- ant iquark Fock state of the longitudinal (L) 
electro- weak boson 2 . For Cabibbo-favored transitions the Fock state expansion reads 

\W£) = ty cs \cs) + ^ ud \ud) + (5) 

where only ud- and cs-states (both vector and axial-vector) are retained. 

In the current conserving eDIS the Fock state expansion of the longitudinal photon con- 
tains only S-wave qq states and ^ vanishes as Q 2 — > 0, 

r) ~ 26 Xt _- x Qz(l - z) log(l/er). (6) 

Here r is the gg-dipole size and z stands for the Sudakov variable of the quark. 
In z/DIS the CCNC adds to Eq.(6) the S-wave mass term [11, 12] 

~ [(m ± fi) [(1 - z)m ± zfi]] log(l/er) (7) 

and generates the P-wave component of ty(z,r), 

~t(6 xrx e-* 2X t>Q- 1 (m±rir-\ (8) 

where upper sign is for the axial- vector current, lower - for the vector one and ( = 2 A - for the 
vector current and £ = 1 - for the axial-vector one. Clearly seen are the built-in divergences 
of the vector and axial-vector currents d v V v ~ m — /J, and d v A v ~ m + fi. This LCWF 
describes the quark-antiquark state with quark of mass m and helicity A = ±1/2 carrying 
fraction z of the W + light-cone momentum and antiquark having mass /x, helicity A = ±1/2 
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Figure 1: Two components of F L = F£ s + F^ d at xsj = 10~ 4 are shown by solid lines. The 
S-wave and P-wave contributions to F[ s and F^ d are represented by dotted and dashed lines, 
correspondingly. 

and momentum fraction 1 — z. The distribution of dipole sizes is controlled by the attenuation 
parameter 

e 2 = Q 2 z(l - z) + (1 - z)m 2 + zfx 2 
that introduces, in fact, the infrared cut-off, r 2 ~ e~ 2 . 



3 High Q 2 : ^- symmetric cs-states 

In the color dipole representation [3, 4] the longitudinal structure function Fl(x,Q 2 ) in the 
vacuum exchange dominated region of x < 0.01 can be represented in a factorized form 

Fl(x,Q 2 ) = J^- J dzd 2 r\^(z,v)\ 2 a(x,r), (9) 

where g is the weak charge, aw = g 2 /^ and Gf/V% = g 2 /™ 2 */- The light cone density of 
color dipole states |\l/| 2 is the incoherent sum of the vector (V) and the axial- vector (A) terms, 

|^| 2 = \v\ 2 + \A\ 2 . (10) 

The Eq. (6) makes it obvious that for large enough virtualities of the probe, Q 2 ^> ml, 
the S'-wave components of both F'£ d and P£ s in expansion (4) are dominated by the "non- 

2 for an alternative description of the vDlS structure functions see e.g. [10]. 
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partonic" configurations with z ~ 1/2 with characteristic dipole sizes [13] 

r 2 ~ Q- 2 . 

In the CD approach the BFKL-log(l/x) evolution [6] of a(x,r) in Eq.(9) is described by the 
CD BFKL equation of Ref.[14]. For qualitative estimates it suffices to use the DLLA (also 
known as DGLAP approximation [15, 16]) Then, for small dipoles [17] 

7r 2 r 2 

a (x, r ) ~ —a s (r 2 )G(x,r 2 ), (11) 

and from Eq.(9) it follows that 

FZ d *F c L s *l-as(Q 2 )G(x,Q 2 ), (12) 

where G(x,k 2 ) = xg(x,k 2 ) is the gluon structure function and asik 2 ) = 4n/ j3o\og(k 2 / A 2 ) 
with fa = 11 - 2N f /3. 

The rhs of (12) is quite similar to of eDIS [15, 18] (see [17] for discussion of corrections 
to DLLA-relationships between the gluon density G and F^). Two S-wave terms in the 
expansion (4) that mimics the expansion (5) evaluated within the CD BFKL model of Ref. [19] 
are shown by dotted curves in Fig. I. The full scale BFKL evolution of the uN structure 
function F L (x, Q 2 ) with boundary condition at x = 0.03 is shown in Fig. 2 of Ref. [20]. 

4 Moderate Q 2 : asymmetric cs-states and P-wave dom- 
inance 

The S-wave term dominates F L at high Q 2 3> m 2 . At moderate Q 2 < ml the P-wave 
component takes over (see Fig.l). To evaluate it we turn to Eq. (9). For m 2 c m 2 in 
Eq.(10), 

1^ H^l 2 °c (gf) e 2 K 2 (er) 

, where K\(x) is the modified Bessel function and one can integrate in (9) over r 2 to see that 
the ^-distribution, dF£ s /dz, develops the parton model peaks at z — > and z — > 1 [7]. To 
clarify the issue of relevant dipole sizes we integrate in (9) first over z near the endpoint z — 1. 
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For r 2 from the region 

(ml + Q 2 )- 1 <r 2 <m; 2 

this yields [8] 

This is the r-distribution for cs-dipoles with c-quark carrying a fraction z ~ 1 of the Vr + 's 
light-cone momentum. Thus, the singularity ~ r~ 4 in Eq.(13) together with the factorization 
relation (9) and a(r) ~ r 2 give rise to nested logarithmic integrals over dipole sizes. Indeed, 
in the Born approximation the gluon density G in Eq.(ll) is 

G(x,r- 2 ) « C F N c L(r- 2 ), (14) 

where 

L(P) = -ilog^. (15) 

Notice, that perturbative gluons do not propagate to large distances and \i G in Eq.(15) stands 
for the inverse Debye screening radius, fie — l/R c - The lattice QCD data suggest R c ~ 0.3 
fm [21]. Because R c is small compared to the typical range of strong interactions, the dipole 
cross section evaluated with the decoupling of soft gluons, k 2 < [i 2 G) would underestimate 
the interaction strength for large color dipoles. In Ref.[22, 19, 23] this missing strength was 
modeled by a non-perturbative, soft correction a npt (r) to the dipole cross section a(r) = 
a ptij') + &npt( r )- Here we concentrate on the perturbative component, <J pt {r), represented by 
Eqs.(ll) and (14). 

Then, for the charmed-strange P-wave component of Fl with fast c-quark {z — > 1) one 

gets 

There is also a contribution to F£ s from the region < r 2 < (m 2 + Q 2 )^ 1 

F? - N i\ m Tl Q ,f s{rn 2 c + Q 2 )L(m 2 c + Q 2 ) (17) 
which is, however one L short. Thus, the CD analysis reveals the ordering of dipole sizes 
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Figure 2: The nucleon structure function F 2 at smallest available x B j as measured in uFe 
CC DIS by the CCFR [24] (circles) and CDHSW Collaboration [25] (squares, x Bj = 0.015). 
Triangles are the CHORUS Collaboration measurements [26] of F 2 in vPb CC DIS. Solid 
curves show the vacuum exchange contribution to F 2 . Also shown are the charm-strange 
(dashed curves) and light flavor (dotted curves) components of F 2) dashed-dotted curves for 
the valence contribution to F 2 . 

(m c 2 + g 2 )- 1 <r 2 «m; 2 (18) 

typical of the DGLAP approximation. The rise of F^[x, Q 2 ) towards small x is generated by 
interactions of the higher Fock states, cs + gluons. The DLLA ordering of Sudakov variables 
and dipole sizes in the n-gluon state \csg\g 2 ...g n ) 

X < Z n < ... < Z! < z < 1 (19) 

(m. + Q 2 )- 1 <r 2 <p 2 < ... <p 2 </i G 2 (20) 
results in the density |<I> n+ i| 2 of multi-gluon states in the color dipole space [3] 
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$ n+1 2 = tt(z,r 2 ^ ^ j 

C F a s {Pi 2 ) 1 Pi CVas(Pn 2 i) 1 p 2 _! 

tH Z 2 P2 7T 2 n p^ 

By virtue of (19,20) the csgig2---g n -st&te interacts like color singlet octet-octet state with the 
cross section (CU/Cj^o^pn). Then, making an explicite use of Eqs. (11,14) and (13) we arrive 
at the P-wave component of F L that rises rapidly to small x, 



4 {ml + Q 2 ) L{ml + 02)V ~ % ( 2 V^ m * + g2) ) • (22) 



iV c C F m 2 



In Eq.(22), which is the DGLAP-counterpart of Eq.(3), 



h, 2 {z) ~ exp(2)/v / 27T^ 

is the Bessel function, 

£0r,£; 2 ) = ?pL(A; 2 ) 

is the DGLAP expansion parameter with rj = Ca \og(x /x). 

Additional contribution to P£ s comes from the P-wave cs-dipoles with "slow" c-quark, 
z — > 0. For low Q 2 <C m 2 this contribution is rather small, 



c 

If, however, Q 2 is large enough, Q 2 ^> m 2 , corresponding distribution of dipole sizes valid for 

{ml + Q 2 )- 1 <r 2 <€.m~ 2 

is 

/^, r)l ^|J_. (24) 

The DLLA summation over the s-channel multi-gluon states, results in [5] 

FT - ~^ L(Q 2 ) V -% (2v^^j) . (25) 

Therefore, at high Q 2 3> m 2 both kinematical domains z — > 1 and (Eqs. (22) and (25), 

respectively) contribute (within the DLLA accuracy) equally to P£ s . 



5 Low Q 2 : light quark dipoles and Adler's theorem. 

The P-wave component of F£ d is small because of small factor m 2 /Q 2 , where m q is the con- 
stituent u, d-quark mass. Here we deal with constituent quarks in the spirit of Weinberg [27]. 
This suppression factor, m 2 /Q 2 , comes from the light-cone wave function ^ ud ~ m^Qr)^ 1 
and is of purely perturbative nature. 

In [20] we checked accuracy of the color dipole description of Fl(x,Q 2 ) in the non- 
perturbative domain of low Q 2 making use of Adler's theorem [2], 

F«(x, 0) = &a w , (26) 

7T 

In (26) f n is the pion decay constant, a n is the on-shell pion-nucleon total cross section. 

Invoking the CD factorization, which is valid for soft as well as for hard diffractive in- 
teractions, we evaluated first the vacuum exchange contribution to both a n and F L (x, 0). 
The parameter /„. in Eq.(26) was evaluated within the CD LCWF technique [28, 29]. The 
approach successfully passed the consistency test: TrFl d (x,0)/(f 2 a 7T ) « 1 to within 10%. The 
cross section a w was found to be in agreement with data. However, the value of f n appeared 
to be underestimated. It was found that for m q = 150 MeV, commonly used now in CD 
models successfully tested against DIS data, our F L at Q 2 — > undershoots the empirical 
value of f%a w /ir by about 40% [20], not quite bad for the model evaluation of non-perturbative 
parameters. One can think of improving accuracy at higher Q 2 ~ ml which we are interested 
in. 

Notice, that Adler's theorem allows only a slow rise of F£ d (x, 0) to small x, 

/l\ A soft 

F£ d (x, 0) oc y , (27) 

much slower than the rise of F[ s following from our DLLA estimates. The value of the so- 
called soft pomeron intercept A so f t ~ 0.08 comes from the Regge parameterization of the 
total irN cross section [30]. 
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6 Comparison with experimental data. 

We evaluate nuclear (uA) and nucleon (vN) structure functions within the color dipole BFKL 
approach [19] (for alternative approaches to nuclear shadowing in neutrino DIS see [31, 32, 
33, 34, 35]). 

The structure function F 2 for the vFe and vPb interactions are shown in Fig. 2. From 
comparison with experimental data [24], [25] and [26] we conclude that the excitation of charm 
contributes significantly to F 2 at x < 0.01 and dominates F 2 at x < 0.001 and Q 2 <m 2 . 

For comparison with data taken at moderately small- a; the valence component, F 2va i, 
of the structure function F 2 should be taken into account. We resort to the parameteriza- 
tion of F 2va i(x, Q 2 ) suggested in [36]. This parameterization gives F 2va i(x, Q 2 ) vanishing as 
Q 2 — > which is not quite satisfactory from the point of view of PCAC. The latter requires 
F 2val (x,0)=F 2 p v c a f c (x,0) with 

71 

Here x = m 2 a /W 2 and a^(W) stands for the secondary reggeon contribution to the total 
pion-nucleon cross section that diminishes at high cms collision energy as a p (W) ~ (W 2 ) aR ~ l , 
where ocr ~ 0.5. However, at smallest values of Q 2 ~ 0.2 — 0.3 GeV 2 accessible experimentally 
F 2va i(x,Q 2 ) 3> F 2 p ^f c (x,0), remind, the characteristic mass scale in the axial channel is 
m a ~ 1 GeV. Therefore, the accuracy of F 2va i(x,Q 2 ) of Ref.[36] is quite sufficient for our 
purposes. In Fig. 2 the valence contributions to F 2 are shown by dash-dotted curves. The 
agreement with data is quite reasonable. 

One more remark is in order, the perturbative light-cone density of ud states, \^> ud \ 2 ~ r~ 2 , 
apparently overestimates the role of short distances at low Q 2 (see Ch. 5) and gives the value 
of Ff d (x ) 0) which is smaller than the value dictated by Adler's theorem [20]. This also may 
lead to underestimation of F 2 in the region of moderately small x > 0.01 dominated by the 
light quark current. 
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7 Summary 

Summarizing, it is shown that at small x and moderate virtualities of the probe, Q 2 ~ m 2 , 
the higher twist corrections brought about by the non- conservation of the charmed-strange 
current dramatically change the longitudinal structure function, Fl. The effect survives the 
limit Q 2 — > and seems to be interesting from a point of view of feasible tests of Adler's 
theorem [2] and the PCAC hypothesis. 
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